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The widths of charge balance function in high energy hadron-hadron and relativistic heavy ion 
collisions are studied using the Monte Carlo generators PYTHIA and AMPT, respectively. The 
narrowing of balance function as the increase of multiplicity is found in both cases. The mean 
parton-freeze-out time of a heavy-ion-collision event is used as the characteristic hadronization time 
of the event. It turns out that for a fixed multiplicity interval the width of balance function is 
consistent with being independent of hadronization time. 

PACS numbers: 13.85.Hd, 25.75.-q, 24.10.Lx 



I. INTRODUCTION 

The relativistic heavy ion collision experiments at 
CERN-SPS and especially at the relativistic heavy ion 
collider RHIC in Brookhaven National Lab provide clear 
evidence for the production of a dense matter in the col- 
lision processes Jj . The central question is whether this 
matter is purely hadronic or has been going through a 
quark-parton phase. There exist experimental evidences 
in favor of the existence of a quark-parton phase at the 
early stage of collision processes "5] , but in view of the im- 
portance of this question, further confirmation is needed. 

Recently, the rapidity correlation between oppositely 
charged particles, which has been used in e+e~ ^ 
and hadron-hadron collisions Q to study the hadroniza- 
tion in these processes, is proposed [3 as a measure 
of the hadronization time in relativistic heavy ion col- 
lisions. It is argued that if the system produced in heavy 
ion collisions has undergone a quark-parton phase, the 
hadronization will occur at a later time and, therefore, 
the temperature will be lower and the diffusive inter- 
action with other particles will be lesser than those in 
the direct hadronization without going through a quark- 
parton phase. These will result in a narrower charge 
balance function for a system with quark-parton phase 
than that without such a phase. 

Two heavy ion experiments 0, 01 have measured the 
balance function at various centralities and for differ- 
ent colliding nuclei. A narrowing of the balance func- 
tion is indeed observed with increasing centrality of the 
collision and with increasing size of the colliding nuclei. 
These observations are consistent with the assumption 
that the narrowing of balance function is correlated with 
late hadronization. 

On the other hand, recently it is reported Q that in 
hadron-hadron collisions at y/s — 22 GeV the balance 
function also becomes narrower as the increasing of mul- 
tiplicity. Therefore, whether the observed narrowing of 
balance function in relativistic heavy ion collisions is due 



to late hadronization or is simply due to the multiplicity 
effect is an open question. 

In this letter this question is examined using the Monte 
Carlo generators PYTHIA [l^ and AMPT [ill- The 
former is a standard Monte Carlo generator with string 
fragmentation as hadronization scheme. There is not any 
quark-parton phase in this model and the hadronization 
is almost instantaneous. On the other hand, the latter is 
a "multi-phase" model, with a transport of quark-parton 
before hadronization. 

The results from PYTHIA will be presented in section 
II. The hadronization time in AMPT model is discussed 
in section III, and its connection with the width of bal- 
ance function is presented in section IV. Section V is 
summary and discussion. 



II. THE WIDTH OF BALANCE FUNCTION IN 
PYTHIA MODEL 



The balance function is defined as 
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where, n-| — (Sy), n — \-{Sy) and n++{Sy), n--[5y) are the 
numbers of pairs of opposite- and like-charged particles 
satisfying the criteria that they fall into the rapidity win- 
dow and that their relative rapidity equals 5y; and 
ri_ are the numbers of positively and negatively charged 
particles, respectively, in the interval Yw- 

The balance function _B((5y|Yw) represents the proba- 
bility that the balancing charges are separated hy Sy 0] . 
The mean of Sy 
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is defined as the width of balance function. 

Proton-proton collision events are generated at four 

cm. energies 22, 64, 130 and 200 GeV using 

PYTHIA5.720 generator. The event number for each 
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FIG. 1: The width of fuU-phase-space balance func- 
tion for different multiplicities in p-p collisions at ^/s — 
22, 64, 130, 200 GeV. 



sample is 100,000. The widths (Sy) oo of balance func- 
tion in the full phase space are calculated for different 
(charged) multiplicity bins and plotted in Fig. 1. 

It can be seen from the figure that in this model even 
for p-p collision, where no quark-parton phase is expected 
and the hadronization is almost instantaneous, the width 
of balance function decreases with the increase of multi- 
plicity, i. e. the width of balance function is narrower for 
higher multiplicity. This effect has nothing to do with 
hadronization time. 

On the other hand, by definition balance function mea- 
sures the correlation length between oppositely charged 
particles. For comparison we have calculated the stan- 
dard 2-particle correlation function of oppositely 
charged particles 
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for different multiplicities in p-p collision at cm. energy 
Vs = 200 GeV, for yi = 0, 7/2 = V- The results plotted 
in Fig. 2 show that the width of R is consistent with be- 
ing independent of multiplicity. A possible explanation 
of the width of R is cluster decay. Comparing with the 
definition of balance function, Eq. (1), we see that it is 
the difference between the correlations of opposite- and 
like-charged particles that shows a clear multiplicity de- 
pendence, which is unrelated with cluster decay, and is 
mainly due to the string fragmentation mechanism im- 
plemented in PYTHIA model. 

It can also be seen from Fig. 1 that the width of bal- 
ance function depends on collision energy. For the same 
multiplicity, the higher the collision energy is, the wider 
the width of balance function. 

However, it should be noticed that the full rapidity 
region is wider for higher energy, c/. Fig 3. In order to 
get rid of the influence of the width of rapidity region 
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FIG. 2: The 2-particle correlation function i?"' (0, y) of 
oppositely charged particles as function of y for differ- 
ent multiplicities in p-p collision at cm. energy y/s = 
200 GeV. 
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FIG. 3: The rapidity distribution of p-p collision at ^/s ■ 
22, 64, 130, 200 GeV. 
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FIG. 4: The width of balance function in the rapidity 
region [—3, 3] for different multiplicities in p-p collision 
at = 22, 64, 130, 200 GeV. 
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we calculate the width of balance function in the region 
— 3 < y < 3 for all four energies. The results, presented 
in Fig. 4, show that, when the (average) rapidity density 
is the same, the width of balance function is almost 

independent of energy, especially for high That is, 
in hadron-hadron collisions the width of balance function 
depends essentially only on multiplicity and is consistent 
with being independent on energy. 

How does the width of balance function behave in 
nucleus-nucleus collisions will be studied in the next two 
sections. 
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III. HADRONIZATION TIME IN AMPT MODEL 

The Monte Carlo generator AMPTl.ll is a multi- 
phase transport model, which contains a quark-parton 
transport phase before hadronization. The initial spa- 
tial and momentum distributions of hard partons and 
soft string excitations are obtained from the HIJING 
model. The parton cascade follows Zhang's parton- 
cascade (ZPC) model which is based on two-body 
pQCD scattering with screening masses. When interac- 
tion ceases, the partons are recombined with their par- 
ent strings to form hadrons according to LUND string 
fragmentation mechanism 10] . Then the scatterings 
among the resulting hadrons are described by a relativis- 
tic transport (ART) model 15] which includes baryon- 
baryon, baryon-meson and meson-meson elastic and in- 
elastic scatterings. 

In our calculation the parton cross section is chosen to 
be 10 mb. If the colliding nuclei are large and the energy 
is high, then the parton cascade will last long enough to 
make the par ton distribution arrive at an equilibrium dis- 
tribution [l6j. However, no equilibrium thermodynamics 
has been included in the model. In particular, there is 
not any equilibrium phase transition from parton phase 
to hadron phase. Therefore, there is no unique hadroni- 
sation time for the whole system. Each parton has its 
own hadronisation time, or freeze out time tfr. 

In order to study the correlation, if any, between the 
width of balance function and the hadronization time, we 
use the event mean of parton freeze out time 
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as the characteristic hadronization time for an event, 
where nparton is the number of partons in the event, tf^i 
is the freeze out time of the ith parton. 

The AMPTl.ll (default) generator is utihzed to gen- 
erate Au-Au collision events at ^/snn = 200 GeV. Two 
event samples with 250,000 and 770,000 events, respec- 
tively, are generated for two centralities 6 < 7 fm and 
6 > 7 fm. 

In Fig's. 5 are shown the scattering plots of ifj- vs. 
rich in Au-Au collisions at ^/snn — 200 GeV for the 
two different centralities 6 < 7 f m and b > 7 fm. 



FIG. 5: Scattering plots of <tr vs. Tich in Au-Au collision 
at y/SNN — 200 GeV for two different centralities. 
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FIG. 6: The distributions of event-mean parton freeze 
out time for two different centralities in Au-Au colli- 
sion at y/SNN ~ 200 GeV. 



respectively. It can be seen that in central collisions (6 < 
7 fm) is larger than 5 fm, while in peripheral collisions 
{b > 7 fm) tfr is concentrated at th < 5 fm. That is, 
central collision events hadronize later than peripheral 
ones. 

The distributions of event-mean parton freeze out time 
tfr for the two different centralities in Au-Au collision at 
^/SNN = 200 GeV are shown in Fig. 6. 



IV. RELATION BETWEEN THE WIDTH OF 
BALANCE FUNCTION AND HADRONIZATION 
TIME 

In order to study the correlation between the width of 
balance function and the characteristics of single event — 
event- mean parton freeze out time tfr and/or multiplicity 
rich, each centrality sample is divided into sub-samples 
according to the intervals of mean parton freeze-out time 
tfj- and the resulting sub-samples are further divided into 
sub-samples by different multiplicity intervals. 

The widths of balance function in the rapidity region 
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FIG. 7: The widths of balance function in the rapidity 
region Yyv = [~3j3] for different mean parton freeze- 
out time ifr intervals versus multiplicity rich in Au-Au 
collision at ^^snn ~ 200 GeV. 



= [^3, 3] for different intervals of mean parton freeze- 
out time tfj- and multiplicity rich are shown in Fig. 7. The 
abscissa of the figure is the average multiplicity in the 
corresponding multiplicity interval. 

It can be seen from the figure that in all cases the 
width of balance function decreases with the increasing 
of multiplicity, while in the same multiplicity interval, the 
width of balance function is consistent of being constant, 
independent of the hadronization time. Therefore, to use 
the narrowing of balance function in rclativistic heavy 
ion collision as a measure of hadronization time and as a 
signal of QGP is doubtful. 



V. SUMMARY AND DISCUSSION 



It is found using PYTHIA Monte Carlo that the width 
of charge balance function decreases with the increasing 
of multiplicity in high energy hadron-hadron collisions, 
where the hadronization is almost instantaneous. 

The relation between the hadronization time and the 
width of charge balance function in rclativistic heavy ion 
collisions is examined using the default AMPTl.ll Monte 
Carlo generator. The mean parton freeze out time of an 
event is used as the characteristic hadronization time of 
the event. The narrowing of balance function as the in- 
crease of multiplicity is found to exist also for rclativistic 
heavy ion collisions, while for a fixed multiplicity interval 
the width of balance function is consistent with being in- 
dependent of hadronization time. Therefore, to use the 
narrowing of balance function in rclativistic heavy ion 
collisions as a measure of hadronization time and as a 
signal of QGP is doubtful. 

It should be noticed that AMPT model is a multi- 
phase transport model. There is no equilibrium thermo- 
dynamics included in the model. In particular, there is no 
equilibrium phase transition, and consequently no unique 
hadronization time for an event. To use the average of 
parton freeze out time in an event as the characteristic 
hadronization time of the event is a crude approximation. 
Therefore, the observation made in the present work on 
the independence of the width of balance function on 
hadronization time is only a first step. Further investi- 
gation along this line is needed. 
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